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^ , Fitting whole spectra at intermediate spectral resolution (R = 1000 - 3000), to 

^ ' derive physical properties of stellar populations, appears as an optimized alternative 

to methods based on spectrophotometric indices: it uses all the redundant information 
contained in the signal. This paper addresses the validation of the method and it inves- 
tigates the quality of the population models together with the reliability of the fitting 
procedures. Our method compares observed optical spectra with models to derive the 
age, metallicity, and line broadening due to the internal kinematics. It is insensitive to 
the shape of the continuum and the results are consistent with Lick indices but three 
times more precise. We are using two algorithms: STECKMAP, a non-parametric 
regularized program and NBURSTS a parametric non-linear minimization. 

We compare three spectral synthesis models for single stellar populations: Pegase- 
, HR, Galaxev (BC03) and Vazdekis/Miles, and we analyse spectra of Galactic clusters 

whose populations are known from studies of color-magnitude diagrams (CMD) and 
spectroscopy of individual stars. 

We find that: (1) The quality of the models critically depends on the stellar li- 
L>( ■ brary they use, and in particular on its coverage in age, metallicity and surface gravity 

I of the stars. Pegase-HR and Vazdekis/Miles are consistent, while the comparison be- 

■ tween Pegase-HR and BC03 shows some systematics reflecting the limitations of the 

stellar library (STELIB) used to generate the latter models; (2) The two fitting pro- 
grams are consistent; (3) For globular clusters and M67 spectra, the method restitutes 
metallicities in agreement with spectroscopy of stars with a precision of 0.14 dex; (4) 
The spectroscopic ages are very sensitive to the presence of a blue horizontal branch 
(BHB) or of blue stragglers. A BHB morphology results in a young SSP-equivalent 
age. Fitting a free amount of blue stars in addition to the SSP model to mimic the 
BHB improves and stabilizes the fit and restores ages in agreement with CMDs stud- 
ies. This method is potentially able to disentangle age or BHB effects in extragalactic 
clusters. Full spectrum fitting is a reliable method to derive the parameters of a stellar 
population. 
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1 INTRODUCTION 

The goal of this paper is to validate a method for determin- 
ing the characteristics of a stellar population using an entire 
medium resolution spectrum in the rest-frame optical range 
* E-mail: mina.koleva@obs.univ-lyonl.fr (rather than concentrating on specific features). Special at- 
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tention is given to make an optimal use of the information, 
to look for systematics and to compare the different models 
of stellar populations. 

Stellar populations keep a record of most of the evolu- 
tion of the stellar systems, and by studying Hertzprung- 
Russell (HR) diagrams (or colour-magnitude diagrams, 
CMDs) of galaxies, it is possible to constrain the history of 
the star formation and enrichment since their formation (see 
eg. Grobcl 2005). Unfortunately, these star- resolved obser- 
vations are restricted to nearby galaxies (Local Group and 
a little further) because of the shallow magnitude limits and 
crowding effect. At large distances, only the spectral energy 
distributions (SED) and the line-of-sight (LOS) integrated 
spectra are available. 

Many methods were devised to determine the age and 
the metallicity of stellar systems from their integrated light. 
In this paper we distinguish three approaches: 

(i) SED fitting - sensitive to the general shape of the con- 
tinuum, (usually) made in broad wavelength ranges; 

(ii) Spectrophotometric indices - concentrates on the 
strength or equivalent width of specific features, not sen- 
sitive to the fiux calibration; 

(iii) Full spectrum fitting - use all the information by fit- 
ting all the pixels, but independently from the shape of the 
continuum; 



S EP fitting llCid Fernandes et all l2005l : iPanter et all 

I2OO3I : ISolorio et al.l |2005| ) can be made using colours, 
low or medium resolution spectra. It is performed over 
large wavelength ranges and compares the observations 
to rn odels hke P egase.2 (iFioc fc Rocca-Volmerangel [l997l . 
I1999I ). Galaxev (jBruzual fc Charlotll2003l ) or other" These 
models should take into account the various sources of ex- 
tinction (internal and Galactic). SED fittig constrains not 
only the stellar population but it also informs about the 
dust content. 

Using low- and medium resolution spectra, the classical 
approach is to measure spectropho tom etric indices, like Lick 
l|Fabeilll973l : IWorthev et allll994l ) or iRose I (psJ) indices 
and compare them to model predictions. Spectrophotomet- 
ric indices measure the strength of some spectral features 
expected to trace the abundance of particular elements or to 
correlate with the mean age of stellar populations. However, 
primarily because a galaxy spectrum is naturally broadened 
by the internal kinematics, the spectrophotometric indices 
are defined at low resolution (a few A) and are therefore 
blends of many spectral lines. That is why their responses to 
age and m etallicity are quit e complex and require careful cal- 



ibrations (jKorn et al. 



2OO5I). Attempts to d efine indices at a 



higher resolution (see lLe Borgne et al.ll2004l ) require delicate 
velocity-dispersion corrections or must be designed carefully 
to reduce their sensitivity to the velocity dispersion. A good 
example of high-reso l ution index is the H7 index defined in 
IVazdekis fc Arimotol l|l999l ). Another possibility to improve 
the sensitivity of an index to either age or metallicity is 
to use a composite definition, like for the CaT* (calc ium 
triplet corrected for Paschen lines, Cenarro et al. 2001) or 
MgFe (insensitive to [a/Fe] , [Gonzalez 1993; Thomas ct al. 
I2OO3I ) indices. 

Lick indices are widely used for more than a decade and 
the methods continue to improve with introducing the ef- 
fects of non-solar abundances. However, they use only a few 



small passbands out of the full observed spectral range, and 
although the information in the spectrum is probably redun- 
dant, this does not optimize the usage of the collected light. 
Another limitation of the in dices is that they are sen sitive 
to missing data or bad pixels (|Chilingarian et ahllioOTl ): The 
interpolation of spectral bins affected by detector defects or 
cosmic hits may bias the measurement of the index if the 
observations are close to the sampli ng limit. The indices are 
also a ffected by weak emission lines (|Goudfrooii fc EmsellemI 
Il996l ). 

The full spectrum fitting that we propose combines the 
advantages of spectrophotometric indices and SED fitting: 
It uses all the information, it is insensitive to extinction or 
fiux calibrating errors and it is not limited by the physi- 
cal broadening, since the internal kinematics is determined 
simultaneously with the population parameters. /.From the 
point of view of measuring the internal kinematics, it is an 
exte nsion of the classical methods (eg. optimal template fit- 
ting, Ivan der Marellll994l). The method has been presented 



in several places (|Prugniel et al.ll2003l : IOcvirk et al.ll2006al lbl) 

and its reliability assessed on the basis of various simula- 
tions. It is three times more sensitive, but still consistent 
with Lick indices (iChilingarian et al...2007. : .Michielsen et al] 
l2007h . 

Full spectrum fitting is well suited when the resolution 
is comparable with the physical broadening (10 kms~^ < 
(Tins < 200 kms~^ ). At lower resolution (R = 500 - 1000) 
the spectral features are diluted and the sensitivity of the 
method would rely on the strong spectral features and there- 
fore would become equivalent to Lick indices. At even lower 
resolution, only strong features, like the Balmer break pro- 
vides information on the stellar population. In this case, 
SED fitting using also the shape of the continuum is the 
good solution. 

In its intention to be insensitive to the fiux calibration 
and extinction our full spectrum fitting is close to spec- 
trophotometric indices, but it uses all the spectrum. A possi- 
ble disavantage compared to Lick indices is that the method 
may, in principle, be sensitive to the wavelengt h range and 
to spectral resolution. However, iKoleva et al. (|2006i ) have 
shown that this sensitivity is not critical: the precision of 
the results mostly depends on the total S/N (i.e. integrated 
over all the spectral range), and even excluding the Balmer 
lines does not bias the age determination (only precision is 
reduced). 

To check the reliability of a spectroscopically deter- 
mined SSP equivalent age/metallicity, it is wise to test it 
on globular clusters (GC) and to compare the results with 
the results fr om CMDs determinations. It has been shown 
(|Woff et al.H2007l ) that the full spectrum fitting reproduces 
better the results from the CMDs than the other methods. 

This paper is organised as follows: In Sect. 2 we present 
the population models and the two algorithms implement- 
ing the full spectrum fitting. In Sect. 3, we use these al- 
gorithms to compare different synthetic population models 
and in Sect. 4 we invert spectra of Galactic clusters and 
compare the results with age determinations from HR dia- 
grams. Section 5 discusses the prospects for analyzing stellar 
populations and gives conclusions. 



Ages and metallicities from full spectrum fitting. 
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2 POPULATION MODELS AND FITTING 
METHOD 

The goal of this analysis is to retrieve simultaneously the 
characteristics of the internal kinematics and of the stellar 
population using a medium resolution spectrum that inte- 
grates along the line-of-sight the light from a region of a 
galaxy or cluster (hereafter LOS-integrated spectrum). 

The main characteristics of this method are (1) to be 
insensitive to the shape of the continuum and (2) use all the 
information in the spectrum. Its principle is to compare an 
observation to a population model convolved with a line-of- 
sight velocity distribution (LOSVD). The spectrum of the 
population model must have the same line-spread function 
(LSF) as the observation (i.e. the LSF of the observation 
must be injected into the model). 

In this section we describe the main ingredients of the 
method used to analyse stellar populations. First we present 
the population models and the stellar libraries used to con- 
struct the synthetic spectra. Then we describe the two al- 
gorithms used to retrieve the stellar population kinematics 
and characteristics from the spectra. 

2.1 Spectral synthesis models for stellar 
populations 

The population model, i.e. the spectrum corresponding to 
some formation scenario including a given evolution, is 
constructed using a population synthesis code. There ex- 
ists several such codes, and in th is paper we are using 
Pegase-HR (|Le Borgne et all |2004 ). Galaxev (BC03) and 
Vazdekis/Miles (Vazdekis et al. in preparation). These codes 
differ by their choice of the physical ingredients (IMF, evo- 
lutionary tracks, ...), by the libraries of stellar spectra they 
are using and by implementation details (in particular in- 
terpolations). 

2.1.1 Pegase-HR and ELODIE 3.1 

The P egase-HR codcQ, as described in iLe Borgne et al.) 
l|2004 l, allows a choice of different IMF and different 
physical ingredients. The evolutionary tracks are from 
Padova group (Bressan et al. 1993 and companion papers: 
Padova 94). They are identical to those used for Pegase.2 
l|Fioc fc Rocca-Volmera ngc 1999). For the present work we 
use SSPs computed with Salpeter IMF (with lower mass 0.1 

Mq and upper mass 1 20 Mp). 

At variance with iLe Borgne et al.l l|2004 ). we are us- 
mg a new version of the ELODIE hbrarjQ (version 3.1; 
iPrugniel et~al] bOOTl ) at a resolution FWHM = 0.55 A or 
i? = AA/A = 10000 at A = 5500 A. The ELODIE li- 
brary is constructed from the s pectroscopic archive o f Ob- 
servatoire de Haute-Provenc^fl (|Moultaka et al.l |2004 ) and 
the principle of the data-reduction, flux calibration and de- 
ter mirMti£n_jof_th^^ parameters are described 
in IPrugniel fc Soubiraril l|200ll '). The set of 1962 spectra of 
1388 st ars is unchanged in the c urrent release with re- 
spect to lPrugniel fc SoubiranI (|2004l ) . A comparison between 

^ http:/ /www2. iap.fr/pegase/pegasehr/ 

^ http: / /www. obs.u-bordeauxl.fr/m2a/soubiran/elodie_library.html 
^ http://atlas.obs-hp.fr/elodie 
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Figure 1. Direct comparison between Pegase-HR/ELODIE.3.1 
(in black) and Vazdekis/Miles (in red, thick line) populations of 
10 Gyr, solar metallicity in the region of H & K. The Pegase- 
HR spectrum was convolved by a LSF of 60 kms~^ dispersion 
to match it to the resolution of MILES, the vertical scale was 
adjusted to have the same flux at 4000 A. 



the different theo r etical a nd empirical libraries made by 
iMartins fc Coelhol ()2007b[ ) shows that the coverage of the 
library in effective temperature, surface gravity and metal- 
licity {Teff, logg and [Fe/H]) is one of the best from the 
present libraries and therefore well suited for population 
studies (still, enlarging the library toward cooler stars and 
lower metallicity is desirable, and is currently under prepa- 
ration). 

In the current release, the wavelength range, 3900 — 
6800 A has been extended in the blue and includes the H & 
K lines (Fig.[l|. 

In addition, the data processing and the determination 
of the atmospheric parameters have been also improved in 
several ways. In particular, it was found that the diffuse light 
was under-subtracted by the standard reduction pipeline. 
This light is normally modeled using the inter-order regions, 
but at the edge of the CCD where the first blue orders are 
located the diffuse light displays steep variations and is not 
properly fitted by the model. The under-subtraction of the 
diffuse light affects the blue region of the spectrum, below 
about 4400 A. We applied an ad'hoc correction of this effect 
in order to avoid the full re-processin g from the initial CCD 
images (see lErspamer fc North|[20o3 for a fine tuned full re- 
processing). The blue range, up to about 4400 A is affected 
by this small correction. 

The comparison of a Pegase-HR SSP with 
Vazdekis/Miles, Fig. [T] shows that the models are vi- 
sually consistent even in the blue extremity where ELODIE 
spectra have lower S/N. It also reveals that the flux calibra- 
tion of ELODIE is losing precision below 3970 A. Though 
it is not critical for our full spectrum fltting method, it may 
be of importance for other applications. This degradation 
is due to the sensitivity drop in the blue and to edge effect 
in the calibration process. 

The isochrones used in Pegase-HR are scaled-solar at 
different values of the total metallicity Z and this has two 
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consequences on the modeling and analysis of stellar popu- 
lations: 

(i) It is well established and understood that the metal- 
lic mixture of the Sun is not universal. For example, non- 
solar abundances of the a-elements are surely the explana- 
tion of the miss-fit of a strong spectral feature such as the 
Mg triplet and MgH ban d near 5170 A in elliptical galax- 
ies (|Worthev et al.il 19921 ) . The globular clusters are other 
important environment where the a-elements are enhanced, 
[Mg/Fe] > 0. Therefore, scaled-solar models are not strictly 
appropriate in those cases, because the location of the 
isochrones in the H R diagram depends on t he detailed chem- 
ical composition (ISalaris fc Weisj Il998l : iKim et al.l |2002| : 
ISalasnich et al.ll2000l ). Age determinations may be biased by 
up to 20 per cent. At low metallicity ([Fe/H] < —0.7 dex) 
the location of the isochrones mostly depends on the total 
metallicity, Z, while at higher metallicity the position of the 
giant branch is rather controlled by [Fe/H] (but a simple 
rescaling is not possible) . This limitation will be solved only 
when evolutionary tracks at non solar [a/Fe] will be imple- 
mented in Pegase-HR; 

(ii) The stellar library is constituted from stars of the 
solar neighborhood, and though detailed measurements of 
[a/Fe] are available for only a small fraction of the stars, 
we can assume that it ha s the typical abundance pattern of 
this population (see e. g. I Wheeler. Sneden fc Truranlll989l : 
[a/Fe] increases when [Fe/H] decreases). In Pegase-HR we 
neglected the change of [a/Fe] and of Z when [Fe/H] de- 
creases (the predicted spectra have the correct Fe content 
but have an excess of a-elements at low metallicity). This 
has to be considered when interpreting the metallicities de- 
termined using these models; 



2.1.2 Galaxev and the STELIB library 

The Galaxev models (jBruzual fc Cha riot' commonly 
called BC03) use a variety of evol utionar y tracks: Padova 94 
(as Peg ase-HR), Padova 2 000 l|Girardi et al.i I2OO0I ) and 
Geneva l|Schaller et al.lll992l ). After comparing these differ- 
ent prescriptions, the authors adopted as 'standard' model 
the Padova 94 set, because Padova 2000 would predict ex- 
agerately old ages for elliptical galaxies (the red giant branch 
is 50 - 200 K warmer in Padova 2000). 

These models are build using the STELIB library 
l|Le Borgne et al.ll2003l ). with a resolution ?a 3 A FWHM 
in the spectral range 3200 - 9500 A. STELIB contains 249 
spectra but only 187 stars have measured metallicity and 
can be used to compute the predicted spectra. The coverage 
in the parameters space is much more limited tha t in the 
other recent libraries fsee iMartins fc Coelhci|[2007bl ). 

2.1.3 Vazdekis and the MILES library 

Vazdckis/Milos modclfl ar e based on th e previous models 
from (Vazdekis 1999 : Vazdekis et al.ll2003h . They use Padova 
2000 isochrones, which as we pointed before has hotter red 
giant branch. They cover ages between 0.1 and 17.5 Gyr and 
metallicities between -1.7 and 0.2 dex. The MILES library 

* |http://www.ucm.es/info/Astrof/miIes/miles.html| 
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Figure 2. Fits of BC03 with NBURSTS (a),Vazdekis/Miles with 
NBURSTS (b) and Vazdekis/Miles with STECKMAP (c) SSPs, 
using Pegase-HR as a reference. Each vector represents an indi- 
vidual inversion. The dots are the location of the nominal char- 
acteristics of the population and the tip of the vectors are the 
values returned by the analysis procedure. 



ISanchez-Blazguez et al.l (J2OO6') used for these models have a 
2.3 A spectral resolution (FWHM) and cover the wavelength 
range from 3525 to 7500 A. This library is believed to be 
better flux calibrated than any other. 



Ages and metallicities from full spectrum fitting. 5 





20 


1 1 1 1 1 1 1 1 1 ; 1 


1 ' ' ' ' 1 ' ' ' ' M 




10 






J 









> 
to 


-10 








-20 




1 , < , , 1 , , < 




90 








80 






70 






c- 


60 




1 , , rTN~T-T-r-^-l-' 



4500 5000 5500 6000 6500 
Wavelength (A) 



Figure 3. Error of velocity and variation on instrumental veloc- 
ity dispersion in BC03/STELIB models as a function of the wave- 
length. The dots represent the measured values and the lines are 
the polynomial fits with degrees of 4 and 2 for velocity and for 
dispersion, respectively. 

2.2 Fitting methods 

2.2.1 Pre-treatment: Line spread function matching 

In order to use SSP models to analyse observed spectra, a 
prior requirement is to match the spectral resolution of the 
models to that of the observations. 

The spectral resolution is coarsely characterized by its 
FWHM, but entering in the details, the instrumental broad- 
ening, or line-spread function (LSF, the spectral equivalent 
of the PSF for images), is not necessarily a gaussian and 
varies with the wavelength and position in the field (for in- 
tegral field or long slit spectrographs). Therefore, matching 
the resolution is not a mere convolution. 

Note that we are considering here only the instrumental 
broadening. Observed spectral features are also broadened 
by the internal kinematics of the galaxies or clusters, and 
this physical broadening is measured by our analysis method. 

The operation proceeds in two steps: First we determine 
the LSF and its variation with the wavelength, and then we 
inject it in the model. To be more precise, we do not need 
to determine the intrinsic, but simply the relative LSF, i.e. 
the relative broadening function between the spectrum to 
analyse and the model. 

In principle, when dealing with original observations, 
the LSF can be determined using some template stars or us- 
ing twilight spectra (i.e. solar spectrum). But in the present 
paper where we analyse models and stars clusters from a 
library, we cannot access such calibrations. Therefore, we 
measure the total broadening (instrumental plus physical). 

To determine the relative LSF we first used the pro- 
gram NBURSTS (see below) to identify the best matching 
Pegase-HR SSP, and then we measured locally the LSF in 
10 overlapping segme nts of the spectru m with the PPXF 
progranfl ( CappcUa ri fc Emsellemll2004 ). 

The ability to recover kinematical information distin- 

^ http:/ /www.strw.leidenuniv.nl/'mcappell/idl/ 



guishes our method from others which either fix the sys- 
temic velocity (vsys) and the velocity dispersion (a) or adopt 
values from independent measurements. In the present case 
where we analyse models and clusters, with no or very low 
physical velocity dispersion, we decreased the measured LSF 
to simulate a physical dispersion of 30 kms~^. 

Finally, the population models are convolved by all the 
local LSFs and interpolated between these convolved seg- 
ments to reproduce the wavelength variation of the real LSF. 
The importance of the LSF injection cannot easily be judged 
from the visual inspection of the residuals between the 'ob- 
servation' and the best-fit model, since the fit without this 
correction appears already excellent. But the effect could 
be seen on the maps: the well in the minimum region 
becomes sharper with the LSF injection and the estimated 
error bars decrease. 

To compute the values of the we are using the follow- 
ing formula: x^ = [{observation — model)^ /errors^]/nfree, 
where n free is the degree of freedom and the errors are sum 
between the errors in the models and the errors in the ob- 
servations (errors^ojjejs -I- error-s^i,^). The errors of the ob- 
servations of globula r clusters are taken from the library of 
ISchiavon et al.|[2005l . For the errors of the models we adopted 
S/N=100 for Vazdekis/Miles and BC03 and 150 for Pegase- 
HR (from visual inspection of residuals from the fits). These 
choices affect only the x^ scale which should therefore be 
considered as relative, as values lower than one on some of 
the maps reported in this paper are demonstrating. Finally 
to produce x^ maps we scan the space of age and metallic- 
ity, fitting only the coefficients of LOSVD and multiplicative 
polynomial and fixing the parameters of the SSPs (t,[Fe/H]). 



2.2.2 NBURSTS 

The program NBURST S (see a preliminary description in 
IChilingarian et ai]|2007l ) performs a parametric non-linear 
fit. 

The model used in this paper is a SSP convolved with 
a gaussian LOSVD and multiplied by a polynomial whose 
degree is chosen to absorb any flux calibration systematics 
or effects of extinction: 

Obs(A) = P„(A) X ( SSP{t, [Fe/H],X)® 

LOSVD{v,ys,a,h3,h4)) (1) 

The free parameters of the minimization are the charac- 
teristics of the SSP - age (t) and metallicity {[Fe/H]); the 
characteristics of the LOSVD {vsys, cr, h3, h4), and the n+1 
coefficients of the multiplicative function P„, a linear com- 
bination of Legendre polynomials. 

NBURSTS makes a Levenberg-Marquart minimization 
and can fit a combination of SSPs or complex stellar popula- 
tion (in this paper we are fitting only SSPs). The population 
models are spline interpolated over a t-[Fe/H] grid of models 
(a linear interpolation would not work because the deriva- 
tives must be continuous). It is written in IDL/GDL starting 
from PPXF and uses the MPFIlEI procedure. The fitting 
algorithm appears stable, because the free parameters are 
independent enough, and no additional regularization was 

Markwardt, |http://cow.physics.wisc.edu/~craigm/idl/idl.html| 
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Figure 4. Direct comparison between SSPs spectra from Vazdekis/Miles (black line) and the SSPs of Pegase-HR/ELODIE.3.1 (red, 
thick line). The ages/metallicities of the SSPs are: (a) 7 Gyr/0.0 dex; (b) 2 Gyr/-0.68 dex; (c) 12.6 Gyr/-1.68 dex. The Pegase-HR 
spectrum was convolved with a LSF of 60 kms~^ dispersion to match it to the resolution of MILES. The spectra are normalized to 1. 
In green (O-C) we show the residuals. 
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Figure 5. maps for the comparison between Vazdekis/Miles 
and Pegase-HR SSPs spectra of Fig. H] ((a) [Fe/H] = 0.0 dex 
t = 7 Gyr; (b) [Fe/H] = -0.68 dex and < = 2 Gyr; (c) [Fe/H] = 
-1.68 dex, t = 12.59 Gyr) inverted witli Pegase-HR/ELODIE.3.1. 
The minimum is shown with a red circle. The blue star correspond 
to the age/metallicity of the input SSP. The contours represent 



required. At low signal-to-noise the high order terms of the 
LOSVD can be neglected (e. g. fixed to 0). 

The degree of the polynomial can be quite large. Exper- 
imenting with fits against Pegase-HR models (wavelength 
ranee 3900 - 6800 A) we found that the results were stable 
for degrees between 12 and 50. With n = 50 and a wave- 
length range of 3000 A, the fit can absorb variations in the 
continuum at the scale of about 60 A, which is comparable to 
Lick indices (linearly interpolating the continuum between 
two side bands). The high order terms of the multiplicative 
polynomial are well decoupled from all the kinematical and 
population parameters. The maps with n = 50 are more 
regular and have a sharper and deeper minimum. When two 
local minima exist, it is easier to identify the absolute one. 
In principle, at an even higher degree the polynomial should 
become degenerated with the velocity dispersion and metal- 
licity, but we did not find such an effect (we did not explore 
degrees higher than 100). The main inconvenient of a high 
degree is the sensitivity to spikes of emission lines, which 
should be preliminary masked or clipped using a fit with a 
low degree polynomial. 

It is important to stress that the inclusion of a 
multiplicative continuum differs from the 'rectification' of 
the flux calibration used in some o ther algorithms, like 
iMathis. Chariot fc BrinchmannI (|2006l ) ■ where the spectrum 
and the model are divided by their top-hat filtered versions 
(width of 500 A). Since both the observation and the model 
are similarly filtered the results should not be biased but 
the sensitivity to the broad molecular bands of MgH, TiO 
or CN is probably reduced. At variance, our approach does 
not affect the real spectral features, it is rather similar to 
the m ethod of Simien & Prugniel (like in lSimien fc Prugniell 
120021 ) that renormalises the observations using a polynomial 
fitted on the ratio between the observation and the best 
model and iterate until convergence. 



2.2.3 STECKMAP 

STECKMApIII (STEUar Content and Kinematics via Maxi- 
mum a Posteriori) is a Bayesian method to interpret galaxy 
spectra in terms of stellar content. It uses a non-parametric 
formalism, i.e. no shape is assumed a priori for the vari- 
ables, for instance the stellar age distribution. Its principles, 
along with extensive simulations are given in lOcvirk et al.l 
(2006al lbl). The algorithm coded in yorick is now reason- 
ably user-friendl y and the package contains some plugged- 
in SSP models (Bruzual & Charlod l2003l : iLe Borgne et all 
2004; Gonzalez Dclgado ot al. 20o3). 

STECKMAP determines age and metallicity distribu- 
tions (actually weights of the SSPs of the basis). In order 
to compare with the results of NBURSTS, we compute the 
luminosity weighted age as: 



(2) 



where ti is the age of the i-th SSP of the basis and Xi is 
the weight (i.e. fiux fraction) of this component. Similarly, 
we derive the luminosity-weighted metallicity as: 



[Fe/H]i, 



(3) 



^ |http;//astro.u-strasbg.fr/Obs/GALAXIES/stecmap_eng.html| 
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Table 1. Compariso n of the ages/metallicities obtained with our method and those from the literature. Col. 1: ID of the spectra in 
ISchiavon et al. 1 ll2005l ') library consisting in (i) name of the cluster; (ii) (a,b,c) for observations in different nights; (iii)(l,2) for extraction 
in different apertures. Col. 2: The ages, metallicities and corresponding errors from full spectrum fitting with Pegase-HR; Col. 3: like Col. 
2 but fitting also a fraction of h ot stars. Col. 4: ages/metallicities/HBR from the literature. The HBR = (B-R)/(B+V+R), is taken from 
the Harris catalog llHarri£lll996l) 



Pegase-HR Pegase-HR+hot stars literature 

Name age [Gyr] [Fe/H] age [Gyr] [Fc/H] age [Gyr] [Fc/H] HBR 
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These luminosity weighted quantities are expected to 
be comparable to the SSP-equivalent values determined by 
other procedures, including NBURSTS, and actually the 
two types of measurements are often confused in the lit- 
erature. Note however that 'luminosity-weighted' and 'SSP- 
equivalent' are not formally the same. 

In the present paper we are interested in comparing 
simple populations and will therefore not consider the de- 
tails of the full stellar age distributions recovered with 
STECKMAP. 



3 ANALYZING ONE POPULATION MODEL 
WITH ANOTHER ONE AS A REFERENCE 

In this section we will compare different models to check 
their consistency. We will invert SSPs produced by a partic- 
ular population synthesis code using a grid of models made 
with a different code, and compare the retrieved population 
parameters with the nominal ones. 

We restricted the comparison to intermediate-age and 
old populations: 600 Myr - 15 Gyr. We do not discuss the 
younger ages because we do not have reference spectra of 
real clusters for an external validation. Note however that we 
noticed a counter-intuitive dependence of the strength of the 
metallic lines in young Pegase-HR models. At a given age, 
lines like Ca at 4227 A become stronger when the metallicity 
decreases, as a result of the decreasing temperature of the 
isochrones. A validation of the analysis of young population 
should be carried on in a separate future work. 

We also excluded the highest bin of super-solar metal- 
licity generated by Pegase-HR which is absent in BC03 and 
Vazdekis/Miles models. 

The metallicity range of the comparisons is —1.7 < 
[Fe/H] < 0.4 for BC03 and -1.7 < [Fe/H] < +0.2 for 
Vazdekis/Miles. For Vazdekis/Miles and BC03 we use the 
SSPs given by the authors. The fit was made starting from 
4000 A to be consistent with the fit of the globular clusters 
(see Sec. 4). 

Figure [2] presents the result of the inversion of both 
BC03 and Vazdekis/Miles SSPs using a grid of Pegase-HR 
models. We performed the fit using the two inversion algo- 
rithms NBURSTS and STECKMAP. 

The BC03 vs Pegase-HR comparison (Fig. [5^) reveals a 
reasonable consistency for solar metallicity, and for a wide 
range of ages (from 900 Myr to 11 Gyr). For either sub- 
solar or super-solar metallicity the inversion by Pegase-HR 
restores characteristics converging toward solar. At metal- 
licity of [Fe/H] = —0.7 dex the inversion is acceptable, but 
degrades again in the lowest metallicity bin. 

The reasons for these discrepancies can be found in 
the details of the computation of the BC03 spectra. Indeed, 
the 187 stars of the BC03 library are strongly concentrated 
around solar metallicity (1/3 of them have —0.2 < [Fe/H] < 
0.1) as in most of the other libraries. But the small number 
of stars remaining at non-solar metallicity makes the mod- 
eling hazardous. 

To generate a SSP spectrum, BC03 uses the STELIB 
spectra which are the closest to the isochrone sampling 
point. In particular, for each metallicity they choose the 
stars having a compatible [Fe/H]. At solar metallicity the 
HR diagram is sampled by 69 stars with close-to-solar metal- 
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Figure 6. Relations for the GCs. (a) - comparison be- 
tween the metalliciti es found with NBURSTS and those from 
ISchiavon et al.l ||2005| ) (spectroscopy of individual stars). The dot- 
ted lines represent 1 a deviation. The green star represents M67; 
(b) - The upper panel shows the histogram of the ages. The dotted 
lines represent the median age for the clusters; The lower panel 
compares the ages and metallicities that we find and the those 
from the literature. The circles are showing the position from the 
literature and the end of the vectors are our estimations. Every- 
where blue is for the clusters with BHB and red is for the other 
clusters. 



licities. But for [Fe/H] = 0.4 dex, the 41 stars that are used 
have a mean metallicity which is only slightly super-solar. 
Also, the 69 stars used for the [Fe/H] = —0.4 dex bin span 
a metallicity range from solar to [Fe/H] = —0.7 dex with 
an average close to [Fe/H] = —0.2 dex. Not surprisingly, the 
inversions of the corresponding SSPs in the age range 0.5 
to 10 Gyr reproduce this convergence toward solar metal- 
licity. The [Fe/H] = —0.7 dex bin is represented by only 22 
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Figure 7. Same as Fig. [6] but the fit was made adding hot stars 
component to the SSPs models to account for the contribution of 
BHB. 



stars which are actually well distributed around this median 
value, and the comparison with Pegase-HR is more consis- 
tent. The last metallicity bin of BC03, [Fe/H] = -1.7 dex 
is computed from 25 stars spanning a large range in metal- 
licity, and the comparison with Pegase-HR follows the ex- 
pected discrepancy. 

A direct comparison between BC03 and Vazdekis/Miles 
gives the same results, hence confirming this analysis. An- 
other weakness of the BC03 models was revealed by our 
analysis. The inaccurate wavelength calibration of STELIB 
results in velocity errors up to 40 kms~^ (shifts of 0.6 A). 
The wavelength dependence of the velocity error can be ex- 
pressed as: Sv{\) = 8248.5 - 6.5433A + 0.001906A^ - 2.419 x 
10"°'^A^ -I- 1.128 X 10~"A"*, and the change of the velocity 
dispersion: ai„s{X) = 346.620 - 0.0855A + 0.000006A^ A is 
in A, Sv and atns are in kms~^.Sv is relative to ELODIE.3.1 



(whose wavelength calibration is very precise), and ains is 
the absolute instrumental velocity dispersion deconvolved 
from the LSF of Pegase-HR (13 kms -^). The LSF for the 
different SSPs from BC03/STELIB was found consistent de- 
spite the inhomogeneous origins of the spectra. 

Contrasting with the previous, the comparison be- 
tween Vazdekis/Miles and Pegase-HR (Fig. [l]3,c) is con- 
siderably more consistent (the MILES library counts al- 
most as many stars as ELODIE.3.1 and has a compara- 
ble coverage of the space of parameters). Pegase-HR is us- 
ing the Padova 94 tracks and Vazdekis/Miles Padova 2000 
whic h, should result in overes timating the SSP-equivalent 
ages (|Bruzual &: Charlotll200j ). The difference between the 
two sets of isochrones is compatible with the deviations be- 
tween the two models on Fig. [2l3,c. 

Figure|3]compares Pegase-HR and Vazdekis/Miles spec- 
tra for three representative regions of Fig. [2Jd,c. The two 
models at the same age and metallicity are over-plotted with 
only adjustment of the velocity dispersion and multiplicative 
polynomial. The residuals are of the order of 1 per cent, even 
for the lowest metallicity bin, where the S/N of the libraries 
is lower because of the few stars in that region. 

We also studied the dependence of the solution on the 
initial guesses given to NBURSTS in order to identify cases 
where the solution is trapped by a local minimum. In gen- 
eral, the solution is stable for a very wide range of initial 
guesses and the maps (Fig. [5)l confirm that there is an 
unique global minimum. 

Only in the region of old age (>10 Gyr) and low metal- 
licity ([Fe/H] ~ -1.7) the topology of the map becomes 
more complex and the sensitivity to the guesses critical. In 
this regime the age dependence of most features is small and 
irregular. It corresponds to the region where the Lick bi- 
indices grids, like H^ vs. MgFe, fold on themselves, making 
age and metallicity determinations ambiguous. This region 
is of importance for the study of globular clusters, and we 
note that unlike Lick indices our method can determine un- 
ambiguously the age if care is taken to avoid the secondary 
minimum (see Sec. 4.1). Nevertheless, we should also stress 
two limitations of the models in that region (i) the stellar li- 
braries are only sparsely populated and (ii) the metallicities 
are poorly sampled by the isochrones, the 2 values provided 
in Padova 94: [Fe/H] = -1.7 dex and [Fe/H] = -0.7 dex are 
considerably spaced, and intermediate models are interpo- 
lated. 

To estimate the effect of the interpolation of the mod- 
eled spectra between the metallicities we computed spectra 
with intermediate metallicities with Pegase-HR (perform- 
ing linear interpolation) and we inverted them with the 
standard Pegase-HR grid (NBURSTS makes spline inter- 
polations). The errors are generally negligible (0.009 dex 
or 2 per cent on age, 0.006 dex on [Fe/H]), except at low 
metallicity, i.e. between the two widely separated sets of 
isochrones: [Fe/H] — —1.2 dex, where the errors are (0.09 dex 
or 20 per cent in age and 0.02 dex on [Fe/H]). 



4 ANALYSIS OF GALACTIC CLUSTERS 

We have shown that (1) full spectrum fitting can recover 
SSP-equivalent ages and metallicities and (2) the two recent 
population models are consistent with each others (BC03 
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has the disadvantage of using sparsely populated library), 
we want now to assess the reliability of this analysis. For 
this purpose, we will analyse high quality spectra of Galactic 
clusters for which studies of deep color-magnitude diagrams 
are separately available. 

First, we will analyse the 60 spectra from the 
ISchiavon etHI l|2005h library of 40 GCs covering a wide 
range of metallicities. Sec ond, we will fit the sp ectrum of 
the open cluster M67 from [Schiavon et al.1 (|2004al '). The lat- 
ter object will complement in the high metallicity and inter- 
mediate age regime the tests made on the globular clusters 
(M67 has solar metallicity and is « 4 Gyr old). 

4.1 Globular clusters 

The sample of GCs from lSchiavon et al.l (|2005l ) was observed 
at a resolution of about 3.1 A. The observations in the range 
3350 - 6430 A were done with a long slit spectrograph in 
drift scan mode in order to integrate the population within 
one core radius. The mean S/N varies from 50 to 240 A"^ 
depending on the cluster. 
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Figure 8. Light fraction of hot stars used by the best-fitting 
model as a function of the horizontal branch ratio of the clusters. 
As expected, the fraction of hot stars increases with the HBR. 



4.1.1 Analysis 

The fitting was performed in the wavelength range 4000 - 
5700 A. We shortened the original spectra to this range be- 
cause (i) the red part of the observations has poorer quality, 
due to the correction of telluric lines and its inclusion does 
not improve the quality of the fit {x^ is not decreased), and 
(ii) when we use also the blue region of the spectrum, in- 
cluding in particular the H & K lines, the quality of the 
fit is decreased and the results are biased by 1 Gyr toward 
younger ages. The S/N in this region is lower in both the 
models and the observations, but this does not explain the 
larger • The H & K lines are deeper in the models than in 
the observations. Figure [1] shows that the Pegase-HR models 
are still consistent with those of Vazdekis/Miles, and the in- 
version of the GC observations with the latter reproduce the 
same problem. If the mismatch comes from the observations, 
we can either suspect a data-reduction problem (subtraction 
of diffuse light) of a characteristics of the cluster population 
(Ca overabundance; but Ca4227 is well fitted). We did not 
investigate further this question and we excluded this blue 
range. 

We used the same grid of Pegase-HR model and multi- 
plicative polynomial as in the previous section. 

The instrumental velocity dispersion measured from the 
calibration lamp by Schiavon et al. (2005) varies with wave- 
length. The instrumental velocity dispersion is quite con- 
stant above 5000 A at about ai„B = 75 — 80 kms~^ and rise 
up to over 100 kms~^ at 4000 A. Since this is a sensitive 
point of our method we re-determined the LSF as explained 
in Sect. 3. We found the wavelength calibration to be pre- 
cise and the high order moments insignificant, so we only 
corrected the variation of (Tins using the median of the 60 
individual LSFs (one per observed spectrum). 

When analysing the spectra with NBURSTS we found 
that at low metallicity the algorithm happens to stop some- 
times at a local minimum. We solved this problem starting 
from different guesses covering ages from 5 to 15 Gyr and 
appropriate metallicities and we chose the deepest (hope- 
fully absolute) minimum. We discarded solutions older than 



15 Gyr, except for NGC 5986 for which there is no minimum 
younger than that. 

At low metallicity the spectral features are not very 
sensitive to age (see Sec. 3, Fig. [5]). In the bi-index grids, this 
corresponds to the region where the different isochrones are 
very close and even cross each-others. In our analysis, this 
results in the local minima found for some clusters. Choosing 
carefully the initial guess (we are making several fits from a 
grid of guesses) , it is possible to determine the real solution. 

The fits for all the observations are presented in graph- 
ical form in Appendix A. The maps in age and metal- 
licity allow to understand the degeneracies and convergence 
issues. At the minimum, the x'^ a-re in the range 0.63 to 2.88 
and the rms residuals in the whole wavelength range are be- 
tween 0.6 per cent and 1.9 per cent, x^ values lower than 
one are most probably due to overestimation of the errors 
in the models. 

A close examination of the residuals reveals that the 
bottom of the modeled Balmer lines are typically 1.9 per cent 
deeper in the best-fitting models than in the observation (in- 
tensity in the center of the line). Forcing the fit to reduce this 
departure (i.e. increasing the weight of the Balmer lines) typ- 
ically leads to an increase of the age by 1 Gyr (i.e. 10 per cent 
or 0.04 dex). The origin of this mismatch may be either in 
the models or in the observations. We have the same residu- 
als when fitting with the Vazdekis/Miles models, therefore, 
if it is not an observational artifact, the reason would have 
to be investigated in the physical ingredients (evolutionary 
tracks and IMF). Otherwise, an under-subtraction of the 
background, could explain the residuals. 

Two other features are neither fitted precisely: The Mg 
triplet and MgH band around 5175 A and the CN band 
near 4160 A. This is clearly an effect of the difference be- 
tween the abundance pattern of the library and the one of 
the clusters which are a-enhanced compared to the Sun. 
At high metallicity, the library is scaled-solar and therefore 
clusters like NGC6528, NGC6553 or NGC6624 have deeper 
Mg than the best-fitting model. At low metallicity, the li- 
brary is also a-enhanced (it h as the abundance pattern of 
the solar neighborhood; see e. g. l Wheeler. Sneden fc TruranI 
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1 19891 ) and the clusters are well fitted (rms ~ 0.1 per cent). 
Figure Al shows the effect of the abundances for different 
metallicities. Prugniel et al. (2007) have shown that using a 
semi-empirical library it is possible to improve the fit and 
determine the enhancement. 

The CN band near 4150 A is often miss-fitted in the 
GCs because non-solar abundances of C and N. Therefore 
we masked the region 4145 — 4155 A which is about half the 
central bandpass of the CNl and CN2 Lick indices. 



4.1. 



Results 



The results of the fits are presented in Table 1. On Fig. [B^ 
we compare the value of [Fe/H] fou nd with our method and 
the one from lSchiavon et all l|2005h : the standard deviation 
between the two series of determination s is 0.17 dex. The 
reference ages are taken mainly from iSalaris fc Weisd |2002| : 
ISantos fc Piattil l2004l : iKrusberg fc Chabove j |2006| . for the 
clusters for which we did not find age reference we adopted 
values of 10 Gyr. 

The distribution of the ages (Fig. [GJj) appears more 
spread than those determined from CMD analyses for which 
the dispersion is ~ 2 Gyr. The effects that can alter the age 
determinations are (1) contamination of the spectrum by 
field stars, and (2) the horizontal branch morphology. 

Contamination by field stars may only be important 
in some rare cases and in particular for the bulge globular 
clusters (like for example NGC 6553). 

The horizontal branch morphology or existence of blue 
stragglers are well kno wn to affect the spectroscopic deter- 
mination of the ages |de Freitas Pacheco fc Barbuvl 1 19951 : 
iLee et al.|[200ol : iMaraston fc Thomas! I2OO0I '). The h"orizontal 
branch morphology, i.e. the colour distribution of the stars of 
the horizontal branch, primarily depends on the metallicity: 
The most metal poor clusters have generally the larger blue 
extension of the HB. But it presents b road variations (the 
so-ca lled second parameter problem, see iRecio-Blanco et al.l 
I2OO6I for a recent dicussion) weakly correlated with other 
properties of the clusters (density, luminosity, ...). This pres- 
ence of old blue stars (Helium burning cor e) mimic an in- 
termed iate / young age population. However. ISchiavon et al.l 
l|2004bl ). revisiting this question with the spectra of the 
database used in the present paper, have shown the abil- 
ity to constrain the BHB contribution from the integrated 
light spectra using Fe sensitive line and the ratio between 
H^ and H^. 

On Fig. [6] and [7] the clusters are separated in two 
groups according to the value of the horiz ontal branch ratio, 
HBR = {B- R)/{B + V + R) < -0.5 (|Harrislll996l ). We 
adopted the limit between the red and the blue HB morphol- 
ogy clusters at HBR > —0 .5 in order t o attach clusters like 
NGC 2808 (HBR = -0.49. lHarrislll996D known to have blue 
tail l|Castellani et al.ll2006l ) to the 'blue' group. Since our 
models do not reproduce the HB morphology, the clusters 
with blue HB tend to have younger ages. 

We included the effect of blue HB as an additional com- 
ponent in our model. We fitted a non-negative linear com- 
bination between the SSP(t, [Fe/H]) and a set of stars in 
the temperature range from 6000 to 20000 K. This new fit 
returns also the contribution of each star in the best-fitting 
mix. These additional parameters are adding more parame- 
ters to the fit. The consequence is that the space of param- 



eters is more degenerated. However, the overall chi2 values 
are a little bit smaller when hot stars are added. Figure [7] 
shows that the ages of the blue clusters are moved by about 
4 Gyr (0.20 dex), while the other clusters are less affected 
(0.07 dex). Only 5 (from 40) blue HB morphology clusters 
are left with an age younger than 8 Gyr. Also the metal- 
licities are in better agreement with the metallicities from 
ISchiavon et al. 1(2005), the standard deviation is 0.14 dex. In 
Fig. [5] we are showing the dependence of the light fraction 
of the hot star in the best-fitting model on HBR. There is 
a tendency to have have a larger fraction of hot stars when 
the HBR increases (e.g. more blue stars). 



4.2 M67 

The Galactic open cluster, M67, is an interesting case to 
complete the comparison with the clusters since it con- 
tains an intermediate age and metal-rich population quite 
similar to many diff use elliptical galaxies (dE; see eg. 
iMichielsen et al.^l2007^ or to M32. 

W e are using the spectrum presented in jSchiavon et al.l 
l2004al 'l and assembled from spectra of individual stars in 
the clusters (and of some nearby dwarfs to replace the M67 
dwarfs which are too faint to be observed). The Pegase-HR 
fit of M67 is shown in Fig. H 

The SSP-equivalent age 2.90 ± 0.02 Gyr and metaUic- 
ity [Fe/H]= -0.033 ± 0.006 dex are close to those deter- 
mined by spectrophotometric indices from the same spec- 
trum (3.5 ± 0.5 Gyr and [Fe/H]= 0.00 dex. ISchiavon et all 
|2004a) and compati ble with the characteristics obtained 
from CMD analysis jSalaris fc Weissll2002l : ISantos fc Piattil 
l2004h . 

In this high-metallicity regime, the models are of the 
best quality as also shown in the comparisons of Sect. 3, 
and the stellar libraries sample ideally the corresponding re- 
gion of atmospheric parameters. For low luminosity elliptical 
galaxies, we can consider our full spectrum fitting method 
as highly reliable. 



5 CONCLUSIONS 

We investigated the consistency and the reliability of the de- 
termination of SSP-equivalent age and metallicity of stellar 
populations using full spectrum fitting. Our method, suited 
to analyse medium resolution {R > 1000) spectra, compares 
an observation to a model to derive the parameters of the 
stellar population and the kinematical shifting and broad- 
ening. 

We found a precise agreement between the Pegase-HR 
and Vazdekis/Miles models. For metallicities greater than - 
1.0 dex and ages between 600 Myr and 15 Gyr the metallici- 
ties of the models agree to 0.064 dex and ages to 15.5 per cent 
(0.060 dex). At low metallicity ([Fe/H]< -1.4 dex) and in- 
termediate to old ages, the age sensitivity is poor, and care 
must be taken to avoid local minima. 

The BC03 models present systematic biases that are 
probably due to the poor metallicity coverage of the STELIB 
library. The main pattern is that they 'overestimate' the 
metallicity of their models for [Fe/H] > dex and they 
'underestimate' the metallicity of their models for —1 < 
[Fe/H < 0. 
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Figure 9. The upper panel displays the spectrum of M67 dSchiavon et al 1 l2004al '). The lowed panel shows the residuals from the fit with 
Pegase-HR and NBURSTS. Note that the vertical scale for the residuals is ten times expanded with respect to the plot of the spectrum. 



Turning to real data, we found a good agreement 
between the SSP-equivalent age and metallicity of M67 
(2.9 Gyr and —0.03 dex) with those derived from spec- 
trophotometric indices and from CMD analysis. Regarding 
the globular clusters, the miss-match of Mg and CN due 
to non-solar abundances does not prevent a reliable deter- 
mination of [Fe/H]: The co mparison with metallic ity mea- 
surements from the paper of ISchiavon et all l|2005h indicate 
rms([Fe/H]) = 0.14 dex. 

The main difficulty in the interpretation of the inte- 
grated spectra of globular clusters, is the determination of 
the ages. In particular the horizontal branch morphology, 
or presence of blue stragglers can significantly alter the 
ages to misleadingly young value. The Pegase-HR SSPs do 
not model the horizontal branch morphology and for blue 
morphology the fits with NBURSTS returns ages under- 
estimated by up to 30-50 per cent. With STECKMAP 
the blue horizontal branch are identified as a young 
300 Myr) burst. For the other clusters, we find a mean age 
of 10.4 Gyr with a dispersion 2.4 Gyr. Including hot stars in 
the NBURSTS fit to represent the blue horizontal branch, 
the age of the blue clusters are restored correctly and the 
red clusters are less affected. This ability to robustly mea- 
sure the age of the clusters is potentially very interesting for 
the study of extra-galactic clusters. 

In a conclusion, as long as the population can be consid- 
ered as a single burst, the full spectrum fitting can reliably 
be used to determine ages and metallicities. It appears par- 
ticularly robust for the intermediate age populations of high 
metallicity (1 Gyr< t < 10 Gyr, -1< [Fe/H] < 0). 

The full spectrum fitting uses the redundancy of the 
information distributed along all the spectrum and it is not 
possible to identify particular features which constrain either 
the age or the metallicity. Surely, the Balmer lines a re the 
best age indicators, but we have shown ,Kolcva ct aj] ()2006l ) 
that the age can be retrieved even if they are masked. In the 
present work each wavelength point is weighted according 
to its uncertainty and we did not intend to modulate the 
weight of the specific spectral feat ures ac cordin g to their 
age or metallicity sensitivity (as in lPanter et al.l i2003|). We 
may in the future investigate this aspect. 

An other important issue for this method is the mis- 
match due to abundance differences, and in particular to 
Q-element enhancement. The present population models 



are tied to the abundance pattern of the solar neighbor- 
hood and cannot fit the high metallicity globular clus- 
ters and the ellipt ical galaxies . But p urely theoretical li- 
brarie s (see Martins fc Coelhd l2007al ). or semi-empirical 



ones (|Prugniel et al.ll2007l ) will soon allow to build models 
with variable abundances. 
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